Theory of Magnetic Fluctuations in Iron Pnictides 
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Magnetic fluctuations in an unconventional superconductor (U-SC) can distinguish between dis- 
tinct proposals for the symmetry of the order parameter. Motivated thereby, we undertake a study 
magnetic fluctuations in Iron pnictides, tracking their evolution from the incoherent normal, pseu- 
dogapped metal, to the U-SC state. Within our proposal of extended-s-plus s xy inplane gap with 
proximity-induced out-of-plane line nodes, (i) we describe the evolution of the spin-lattice relax- 
ation rate, from a non-Korringa form in the normal state, to a power-law form in the U-SC in good 
agreement with experiment, and (ii) we predict a sharp resonance in the U-SC state along (7r,7r), 
but not along (n/2,0), along with modulated c-axis intensity in inelastic neutron scattering work 
as a specific and testable manifestation of our proposal. 

PACS numbers: 74.70.-b, 74.25.Ha, 76.60.-k, 74.20.Rp 
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The precise mechanism of unconventional supercon- 
ductivity (U-SC) in the recently discovered Iron Pnic- 
tides (FePn) is presently a hotly debated issue [l[ . While 
many physical responses are reminiscent of cuprates Q , 
FePn are metals, albeit presumably proximate to a Mott 
insulator. Moreover, relevance of all d orbitals in FePn 
considerably complicates determination of the pair sym- 
metry. 

Study of magnetic fluctuations in an U-SC can help 
unearth the symmetry of the SC order parameter, as 
shown by detailed studies for cuprates [5]. In the FePn, 
NMR studies already reveal normal state pseudogap be- 
havior 0] and U-SC. The spin-lattice relaxation rate, 
Tf , shows marked deviation from the linear- in-T Ko- 
rringa form expected from a Fermi liquid and smoothly 
decreases for T < 2T C , where T c is the SC transition 
temperature. At very low T, one finds T^ 1 ~ T n with 
n 7^ 3, 5, indicating line nodes in the SC gap. How- 
ever, other probes reveal anisotropic, albeit fully gapped, 
structure of the in-plane gap function. Thus, extant data 
imply that, either one has out-of-plane line nodes, as 
in Sr 2 Ru04 : @, or disorder effects in a s±-SC produce 
the observed behavior [jj]. In the extended- s wave idea, 
disorder is argued to lift the nodal structure, again giv- 
ing similar behavior pj. The issue is thus controversial: 
while ARPES data are inconclusive regarding existence of 
nodes on electron-like FS sheets [1] , a penetration depth 
study, at least in SmO\- x FeAsF x , shows smooth angular 
variation of the in-plane gap Q . While inelastic neutron 
scattering (INS) work does reveal a low-energy resonance 
structure in the U-SC state for Q = (it, tt) [lOj], more de- 
tailed map of the INS response in q-space awaits future 
work. To date, we are aware of one study where the dy- 
namical spin susceptibility, (c[, u>) , nas been measured 
for the 122 FePn, showing that the SC gap has in-plane 
smooth angular variation and an out-of-plane cos(fc z c) 



and strong [13|, [14[ coupling limits. In the itinerant ap- 
proach, the magnetic fluctuations have been computed 
within HF-RPA. For s± pairing, a sharp resonance for 
Q = (it, it) in INS is predicted below T c jJ3|. while no 
such feature is arises for s, ex-s, and d wave-pairing, or for 
q 7^ Q. To get the power-law-in-T behavior in NMR and 
^SR in the s± idea, it is necessary to consider (strong) 
disorder effects in a two- or four-band model. Again, 
the situation is controversial. For LaFePO, the pene- 
tration depth, A(T) ~ T 12 [15], while a similar study 
on a wide range of samples of different FePn found a 
seeming univers ality in A(T); this mitigates against the 
disorder effects 16]. //this is true, one must consider 
both the NMR and INS data within a theoretical sce- 
nario with out-of-plane line nodes in the SC gap, since 
in-plane nodes seem to be ruled out by extant tunnelling 
data 
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Extant theoretical works have studied these issues us- 
ing effective model Hamiltonians, both in the weak 12 1 



Recently, based on inputs from the correlated normal 
state electronic structure and rigorous symmetry argu- 
ments, we proposed a specific gap function with sizable 
in-plane angular variation (but no nodes) and inter-band 
proximity induced out-of-plane line nodes [18j |. In con- 
trast to the itinerant picture, our proposal is based on 
a strong correlation view of FePn. Here, we investigate 
the NMR and INS response within such a correlated ap- 
proach, using the full, multiband spectral functions for 
all d orbitals. LDA+DMFT can readily access the in- 
termediate coupling regime relevant for FePn |19j. We 
show how our proposal gives a quantitative account of 
the NMR T-f 1 over the whole T range, and makes spe- 
cific predictions with regard to the observation of the 
low-energy dispersive resonance in the INS intensity be- 
low T c . 

The central quantity of interest is the dynamical spin 
susceptibility, x(q, w) — J2 a b Xa&(q, w), where a, b are all 
(i-orbital indices, and q, uj are the momentum and energy 
transfers in INS. Viewing FePn as strongly correlated sys- 
tems with U = 4.0 cV, V = 2.6 eV and J H = 0.7 eV, 
we construct x(q, u>) in terms of the full LDA+DMFT 
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propagators computed in earlier work [18|, [20j. Very 
good quantitative agreement between LDA+DMFT and 
key experiments in both, the normal and U-SC states, 
has been shown there, lending strong support for our 
choice. The prescription is simple: replace the band 
Green functions used in weak-coupling approaches [l2j 
by their LDA+DMFT counterparts. This ensures that 
the dynamical aspect of strong, local, multi-orbital (MO) 
correlations is included from the outset. 

For a MO-system, after replacing the bare G aa (li,uj) 
with G aa (k, W ) = G^ A+DMFT (k,cu) = [co - 



G aa (k , Co') , ^r^pTTTT i a nd introducing the spin opera- 



tor S a ^(q) = \Y. 



(k + q)< (T .^c a , M , (7 /(k), with 



p = x,y,z, the "bare" dynamical spin susceptibility 
reads 



a,b(q^) = -\°a,aa< ' CT W ^ [ G aa (k + q, W + Cj') 
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Xo 



G bb (k, </) + F ab (-k - q, -u - u')F ba (k, w')]- 

Including the ladder vertex in an infinite summa- 
tion of "ladder" diagrams using RPA, the renormal- 
ized magnetic susceptibility, Xa,b(q, w) = bCo.l. bi^) ~ 
^(q)]~\ where Xo,a,b{w) = E q Xo,a,b(q, w) and J(q) = 
Ji(cos(g :c a)+cos(gj / a))-|- J2COs(g x a)cos((j a a), with J% ~ 

t 2 t' 2 

17 " being the frustrated su- 
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perexchange scales in FePn Q. Using Xo,a,h(e) = 
G / rfe/(e)[l - /(e)]W(e) in the RPA series, the NMR 



relaxation rate, T, 



can be now 



expressed in terms of the full DMFT propagators. 
Here, W(e) = Y. a , b [P*a{e)p bb {e) + p ab (e) p ba (e)] M 
and the p a a(e), Pab(e) are the LDA+DMFT local spec- 
tral functions computed earlier [l^|. Also, G = 
2(^ L )(7e7 n ft) 2 (p-)- Finally, our restriction to the non- 
crossing diagrams in the ladder approximation for x(q, u)) 
is an approximation. It is possible that "non-crossing" 
diagrams need to be included in a full description. How- 
ever, for the underdoped cuprates, it has been shown that 
a renormalized "RPA" summation for x(q, u>) with fully 
renormalized one-particle G CT (fc, u), F a> ^ (T (k, u>) gives ex- 
cellent reconciliation of ARPES and INS data [22]. This 
suggests small vertex corrections: while we cannot prove 
why this should be the case, we argue that the good 
agreement we find below is an a posteriori justification 
for neglecting them in our theory. 

The NMR spin-lattice relaxation rate is a measure of 
the local spin fluctuation rate in both phases. For an 
s-wave SC, the coherence factors give the Hebel-Slichter 
(HS) enhancement as a peak in T^ 1 below T c . When the 
"normal" state is strongly incoherent (large ImX(u; = 
Ep) ^ 0, as in our case), or the SC gap has nodes 3, 



the HS peak is absent. But T, 



e A / kT survives for 



in T characterizes an U-SC with gap nodes [1,0]. In the 
normal state above T c , we set F ab (k, to) = 0. This suffices 
for computing the NMR T^ 1 . More work has to be done 
to compute the INS intensity; we will present deatils in 
a separate work. 



2.5 
H 2 

- 1.5 

- 1 
-I 0.5 



-0.5 
-1 

-1.5 




T « T c in an s-wave SC, while a power-law fall-off 
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FIG. 1: Proposed gap function for the 1111-Iron Pnictides. 
The gap function has nearest- (Ai, with ex-s symmetry) and 
next-nearest (A2, with s xy symmetry) neighbor components. 
With A2/A1 = 0.375, no in-plane gap nodes arise in the gap 
function, in agreement with experiment §]. 



However, qualitative remarks about what we ex- 
pect in the INS response are possible without 
a full analysis. The in-plane part, A ab (k) = 
Ai(cos(k x a)+cos(kya))+A 2 cos(k x a)cos(k y a), of our pro- 
posed gap function is shown in Fig. [1] With electron- 
and hole Fermi sheets well separated as in LDA (or 
LDA+DMFT), no in-plane gap nodes are possible, in 
agreement with a host of measurements 0, H, E3] • In- 
terestingly, this leads to A(k + Q)A(k) < for k 
along (0,0) - (7r,7r) and to A(k + Q)A(k) > for k 
near (±7r/2, 0), (0, ±tt/2). This implies, following earlier- 
work [l^ ]. that INS measurements will show appearance 
of a sharp collective "spin exciton" mode in the U-SC 
state at k = (7T,7r), but none for k = (±7r/2,0). Of 
course, incoherent features coming from DMFT prop- 
agators will introduce damping of this mode, but the 
qualitative feature should survive. Since an out-of-plane 
cos(fc z c) component is induced in the full gap function 
due to interband proximity [18| effect, the INS intensity 
should also reflect this modulation in q z . This last predic- 
tion is a consequence of our form of the full gap function, 
and goes beyond previous work [12]. Such a resonance, 
albeit sizably damped, is indeed seen in INS work on the 
122-FePn [10(. Moreover, the cos( q z c) form has also been 
measured by INS on the 122 FePn [llj, but remains to be 
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checked in the 1111 family. Finally, the in-plane angu- 
lar modulation of the gap function is inferred from /xSR 
work on the S'm-based FePn Q. Thus, rationalization 
with extant INS results readily follows directly from our 
proposal for the gap function. 

Next, we discuss the NMR relaxation rate in the nor- 
mal and U-SC states, making detailed comparison with 
experimental work. In the normal, incoherent metal 
state, the reduction in Tf 1 below 200 K 0, [23J indi- 
cates opening of a spin gap, as in underdoped cuprates. 
While the spin gap in cuprates has been identified with 
short-range magnetic correlations in a quasi-2£>, doped 
quantum antiferromagnet, its origin in the multi-band 
FePn is not settled. We emphasize that this behavior is 
observed in the same regime where PES [H and optical 
data 2J| show incoherent charge dynamics, corroborated 
by a linear-in-T resistivity, a T-dependcnt Hall constant, 
and no Drude peak in optics. All these are compelling 
indicators of a strongly correlated metal. We regard this 
as a justification for using LDA+DMFT. 
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FIG. 2: (Color online) T-dependence of the NMR (TiT) -1 
over the full T range for LaOi- x FeAsF x , with x — (black), 
x — 0.1 (red, dotted), x — 0.2 (blue, dot-dashed) and with 
inclusion of U-SC for x = 0.1 (green, dashed). Notice how 
a doping-dependent spin gap around T* ~ 150 K opens in 
the dop_ed case (x = 0.1, 0.2), in good agreement with exper- 
iment 

In Fig. El we show the NMR (TiT) -1 as a function of 
electron doping for LaO\- x FeAsF Xl with = 0.0, 0.1, 0.2. 
Since we do not consider the q = (it, 0) SDW phase, 
the x = 0.0 curve should only be trusted above T = 
T N = 135 K (shown by the black curve in Fig. [2]). With 
x = 0.1,0.2, however, SDW order is destroyed, and U- 
SC emerges at low T. In this range of x, our results can 
validly be compared to experiment, which we now turn 
to do. 

Quite remarkably, a direct comparison with published 



NMR work [J] reveals good agreement between theory 
and experiment around x = 0.1. The absence of the 
(TiT) -1 = const regime is striking. In particular, both 
experiment and our result show a quasi-linear-in-T (like 
T 0.8-o.9) increase m \/ Tl at "high" T > 200 K (see in- 
set of Fig. [3]) . This resembles the high-T precursor of a 
quantum critical system, and corresponds to the "strange 
metal" regime in the T vs x phase diagrams for this 
system [25j. However, as T is lowered, a smooth drop 
in (TiT) -1 around 150 K marks the onset of the grad- 
ual opening up of a spin gap. Given strong frustration 
(J2/J1 — 0.7) in FePn, it is tempting to link this spin 
gap with strong, short-ranged AF correlations, which are 
expected to survive the doping induced destruction of the 
SDW 26!]. We note that the J x - J 2 H3] model has also 
been used to provide a quantitative fit of INS results for 
the undoped 122 FePn, though it is formally valid in the 
strictly localized regime. This is additional evidence for 
a strong coupling picture, since, in the itinerant picture, 
melting of the SDW should yield a paramagnetic Fermi 
liquid at low T with no spin gap, at variance with ob- 
servations. From our results, we estimate a renormalized 
spin gap scale O(150) K in the 1111 FePn. 
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FIG. 3: (Color online) Low T behavior of the NMR Tf 1 on 
a log-log plot (main panel) and on a normal scale (inset). 
Clear power-law behavior without the Hebel-Slichter coher- 
ence peak, in good agreement with experiment [J, is seen. 

In fact, Tf w tanh(0A2T) over almost the whole 
range from low-(T > 15 K) to high T. While this is not 
particularly illuminating, it shows that the "marginal" 
form, Xiod 00 ) — — i s recovered only at high T, and 

is cut off by the spin gap around 200 — 250 K. This bears a 
peculiar resemblance to underdoped cuprates. However, 
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is seen. 



at very low T, a power-law form ~ T 
This is intriguing, and is fit neither by self-consistent 
renormalization theory 28| , nor by any known local non- 
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FL exponents [29[. It could involve several, frustrated, 
nearly degenerate spin fluctuation channels coming from 
the multiband nature of FePn, but we are unable to quan- 
tify this further. 

However, we can still make a few qualitative remarks 
to get more insight. In the strongly correlated metal, 
with a very small "coherent" component in the DMFT 
spectral functions @, HH , "Mottness" underpins the low- 
energy physics. More precisely, when one is close to a 
correlation-driven Mott insulator, the metallic state has 
small density of quasi-itinerant carriers co-existing with 
effectively local moments |2j ■ These latter arise from inte- 
grating out the high energy Hubbard bands in the DMFT 
spectral function, as argued by Baskaran, Si et al. and 
Wu et al. Given the frustrated hoppings characteristic of 
FePn, the spin degrees of freedom are qualitatively de- 
scribed by an effective frustrated J± — Ji Heisenberg-type 
model. In this model, there is a large window in T, be- 
tween Ts dw a- n d T s [261 ] , where lattice translational sym- 
metry is spontaneously broken but the spin rotational 
(SU(2)) symmetry is not. This naturally leads to gener- 
ation of a spin gap, in agreement with observations. Of 
course, as LDA+DMFT shows, the actual situation in 
FePn is somewhat removed from a strictly localized limit 
where the J% — J2 model would apply. However, in view 
of the Mottness, we believe that it still provides a quali- 
tative understanding of the features derived above in the 
full DMFT calculation. 

At T c , there is no HS peak, as seen in Fig. [3] in our 
work, this arises from strong inelastic scattering in the 



"normal" incoherent state [18|, |2Cj (notice that (u>) en- 



ters the DMFT equation for G a {oj) in the SC state, pro- 
ducing strong damping). At very low T « T c , Tf X (T) 
shows a power-law-in-T dependence: T 1 _1 (T) ~ T n , with 
11 = 2.2 — 2.5, qualitatively consistent with observations 
in the 1111 FePn 0], which show neither a T 3 nor a T 5 
law for T « T c . The two-step variation of below 
T c is also reproduced theoretically. The first "step" from 
T c > T > T c /4 is dominantly governed by the larger 
gap component, while the lower-T variation comes from 
the smaller gap component, as expected from an in-plane 
anisotropic gap, while the power-law variation is ascribed 
to out-of-plane line nodes in such a gap. It is still pos- 
sible that disorder effects (which must be treated in the 
unitary limit [l[) will lift the out-of-plane gap nodes, as 
discussed by Maier et al. [13] and give T^ 1 ~ T 3 be- 
havior |6j; this remains to be checked. In our theory, 
the power law behavior arises from the out-of-plane line 
nodes, induced in the gap by an interband proximity ef- 
fect. Since it does not require disorder effects, our conclu- 
sion should be more "universal" [l6j| . Thus, our results 
show how good agreement with the NMR data is derived 
in the whole T range in terms of our theoretical picture 
of an U-SC with proximity induced line nodes, arising 
from an incoherent normal state at T c . 

In conclusion, we have studied the magnetic fluctua- 



tions in FePn, based on a novel theoretical proposal for 
the symmetry of the SC gap function. In a picture where 
U-SC with out-of-plane gap nodes arises from an inco- 
herent, strongly correlated normal state, we have shown 
how the T-dependence of the NMR relaxation rate can 
be nicely understood over the whole T range, from the 
lowest- to "high" T. Moreover, we have argued how the 
specific form of the gap function allows for concrete pre- 
dictions concerning the observation of the collective reso- 
nant peak in INS measurements. Our study provides fur- 
ther support for the strongly correlated nature of FePn 
above T c , and puts our theoretical proposal of an U-SC 
with out-of-plane gap nodes on a firmer footing. 
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